Micromagnetic Simulation of nanoscale Films with perpendicular Anisotropy 
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A model is studied for the theoretical description of 
nanoscale magnetic films with high perpendicular anisotropy. 
In the model the magnetic film is described in terms of single 
domain magnetic grains with Ising-like behavior, interacting 
via exchange as well as via dipolar forces. Additionally, the 
model contains an energy barrier and a coupling to an exter- 
nal magnetic field. Disorder is taken into account in order 
to describe realistic domain and domain wall structures. The 
influence of a finite temperature as well as the dynamics can 
be modeled by a Monte Carlo simulation. 

Many of the experimental findings can be investigated and 
at least partly understood by the model introduced above. 
For thin films the magnetisation reversal is driven by domain 
wall motion. The results for the field and temperature de- 
pendence of the domain wall velocity suggest that for thin 
films hysteresis can be described as a depinning transition 
of the domain walls rounded by thermal activation for finite 
temperatures. 
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I. INTRODUCTION 



We focus on thin ferromagnetic films with high per- 
pendicular anisotropy like CoPt and MnBi which can 
be used as magneto-optical storage media. In these films 
two different mechanisms can be thought to dominate 
the reversal process: either nucleation or domain wall 
motion M . Which of these mechanisms dominates a re- 
versal process depends on the interplay of the different 
interaction forces between domains with different mag- 
netic orientation. In a recent experiment on Co28-P^72 
alloy films |2||j| a crossover from magnetisation reversal 
dominated by domain growth to a reversal dominated by 
a continuous nucleation of domains was found depend- 
ing on the film thickness which was varied from 100A to 
300A. Correspondingly, characteristic differences for the 
hysteresis loops have been found. Similar results have 
been achieved by simulations of a micromagnetic model 
using energy minimization techniques and Monte 

Carlo methods, respectively It is the goal of this 
paper to investigate the influence of the field and the 
temperature on the domain wall velocity for the case of 
magnetisation reversal dominated by domain wall motion 
by a Monte Carlo simulation of a micromagnetic model. 



II. MICROMAGNETIC MODEL 

Co28-P^72 alloy films have a polycrystalline structure 
with grain diameters of 100-250A. For a theoretical de- 
scription by a micromagnetic model jg] the film is thought 
to consist of cells on a square lattice with a square base of 
size L x L where L — 200A and height h of 100A. Due to 
the high anisotropy of the Co^Ptii alloy film the grains 
are thought to be magnetised perpendicular to the film 
only with a uniform magnetisation M s which is set to the 
experimental value of M s = 365kA/m for the saturation 
magnetisation in these systems ||. The grains interact 
via domain wall energy and dipole interaction. The cou- 
pling of the magnetisation to an external magnetic field 
H is taken into account as well as an energy barrier which 
has to be overcome during the reversal process of a single 
cell. 

From these considerations it follows that the change of 
energy caused by reversal of a cell i with magnetisation 
L 2 hM s <Ti with <Ji = ±1 is: 
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The first term describes the wall energy AE W . The 
sum is over the four next neighbors and for S w we use a 
value of S w = 0.0022J/m 2 which is approximately 50% 
of the Bloch-wall energy Sb for this system 0|. The 
reason for this reduction of the grain interaction energy 
compared to the Bloch-wall energy is, that the crystalline 
structure of the system is interrupted at the grain bound- 
ary and also that due to their irregular shape the grains 
are not connected via their complete surface Lh. 

In the second term describing the dipole coupling AEd 
the sum is over all cells. r,-j is the distance between two 
cells i and j in units of the lattice constant L. For large 
distances it is v(crj,rij) = -£-. For shorter distances a 

ij 

more complicated form which is a better approximation 
for the shape of the cells which we consider can be deter- 
mined numerically and was taken into account. 

The third term describes the coupling AEh to an ex- 
ternal field H. 

Additionally, an energy barrier 6{ must be considered 
which describes the fact that a certain energy is needed 
to reverse an isolated cell by domain wall motion through 
the grain (see also 0). We assume that during the re- 
versal process the energy barrier has its maximum value 
LhSb when the domain wall is in the center of the cell, i.e. 
when half of the cell is already reversed. Consequently, 
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the energy barrier which is relevant for the reversal pro- 
cess is reduced to S = max(0, LhSb — ^\(AE W + AE d + 
AEh)\)- The simulations are in good qualitative agree- 
ment with experiments using Sb — 0.0007J/m 2 . 

In order to simulate the CoPt films realistically disor- 
der has to be considered. Obviously, the grain sizes are 
randomly distributed . In the model above this corre- 
sponds to a random distribution of L which can hardly 
be simulated exactly since it modulates the normalized 
cell distance r,j of the dipole interaction. Therefore, as a 
simplified ansatz to simulate the influence of disorder we 
randomly distribute L in the energy term that describes 
the coupling to the external field. Here a random fluctu- 
ation of L is most relevant, since this term is the only one 
scaling quadratic ly with L. In the simulations we use a 
distribution which is Gaussian with width A. Note, that 
through this kind of disorder our model is mapped on a 
random-field model. 

The simulation of the model above was done as in an 
earlier publications ffllj via Monte Carlo methods Q us- 
ing the Metropolis algorithm with an additional energy 
barrier. Since the algorithm satisfies detailed balance and 
Glauber dynamics it allows the investigation of thermal 
properties as well as the investigation of the dynamics of 
the system. For temperatures T — > the Monte Carlo 
algorithm passes into a simple energy minimization al- 
gorithm with single spin flip dynamics, so that also the 
case of zero temperature can be investigated. 

The size of the lattice was typically 150 x 150. The 
dipole interaction was taken into account without any 
cut-off or mean field approximation. 

III. HYSTERESIS AND DYNAMICS 

Fig. [j] shows a simulated hysteresis loop. The loop is 
nearly rectangular. Here, the reversal is dominated by 
domain wall motion Q. Once a nucleus begins to grow 
the domain wall motion does not stop until the magneti- 
sation has completely changed. The nucleation field for 
the simulation is too high compared to the correspond- 
ing experimental results |||| . There are several possible 
reasons for this effect. First, the nucleation field which 
in the case of domain wall motion dominated reversal is 
the field where domain wall motion starts depends on the 
size and on the shape of the nucleus. In an experimen- 
tal situation a nucleus can be very large e. g. a scratch 
while there is no such artificial nucleus in our simulation 
except of the disorder which leads only to nuclei of very 
small size. However, it is not the aim of our simulations 
to calculate the nucleation field accurately. Rather, the 
simulations are thought to contribute to a better under- 
standing of the fundamental properties of the system. 

Since the reversal is driven by domain wall motion, the 
velocity of the domain wall is a central quantity which we 
will investigate in the following. For the determination of 
the domain wall velocity within the simulation we start 
with a system that has a nucleus of circular shape with 
a radius of 19 cells in the center of the 150 x 150 system. 
When we switch on the driving field, from the nucleus 
a domain starts to grow. For the better observation of 



the domain growth, in our flip-algorithm we do not con- 
sider cells that are not connected to the growing domain, 
i. e. we exclude the possibility of additional nucleation. 
Otherwise we have - at least for finite temperatures - the 
problem that spontaneously new nuclei are build by ther- 
mal activation which with increasing radius overlap with 
the original domain. Fig. ^| shows a domain during the re- 
versal. The black region is the reversed domain following 
the magnetic field. It has a circular shape with a rough 
domain wall. From the domain configurations the mean 
radius r of the domains can be determined through the 
area F of the reversed domain as r = \J F/tt, assuming 
that the domain has a circular shape. 

In order to get a deeper understanding of the influence 
of the temperature on the dynamics we simulated the 
reversal for temperatures T = 0, 300, and 600K. Fig. || 
shows - as an example - the r(t) behavior from the sim- 
ulations for T — OK and different fields. For r > 75 
the domain reaches the boundary of the system and - 
consequently - r(t) saturates. For zero temperature the 
domain wall is pinned for lower fields, i. e. after a short 
period of rearrangement of the domain wall the domain 
wall movement stops and the radius remains constant. 
The pinning of the domain wall is due to energy barri- 
ers which follow from the disorder, the dipole field, and 
the intrinsic energy barrier of the single cell. For finite 
temperatures, the domain wall velocity is always finite. 

For 20 < r < 75 the slope of the r(t) curve is approx- 
imately constant and v can be determined by fitting to 
a straight line. Fig. [| shows the dependence of the do- 
main wall velocity on the driving field for T = 0, 300, 
and 600K. For zero temperature there is a sharp depin- 
ning transition Q at a critical field H c from a pinned 
phase with v = to a phase with finite domain wall ve- 
locity. This transition can be interpreted in terms of a 
phase transition with v ~ (H — H c ) 9 for H > H c where 
in our case the critical exponent is 9 m 1, a value which 
is the mean field result for a moving elastic interface |Tc| ] 
in a random field. Also, this value has been observed in 
simulations of a soft spin model with random-fields [ p"l| . 
For finite temperatures the transition is smeared since 
for finite temperatures there is even for H < H c for each 
energy barrier a finite probability that the barrier can 
be overcome by thermal fluctuations. The corresponding 
waiting time can be expected to be exponentially large 
so that for H < H c the domain wall velocity should de- 
crease like Iilv ~ (H — H c )/T. To illustrate this in Fig.|| 
we show the corresponding semi-logarithmic scaling plot. 
As we expect, the data for the two different finite tem- 
peratures collapse for H < H c on a straight line. For 
H > H c thermal activation is obviously less relevant. 
Here, the dynamics is dominated by the zero tempera- 
ture depinning transition. 

To conclude, the results for the domain wall velocity 
suggest that for zero temperature the hysteresis driven 
by domain wall motion can be understood as a depin- 
ning transition of the domain walls. For finite tempera- 
ture the transition is rounded and for fields smaller than 
the depinning field the domain wall movement is domi- 
nated by thermal activation. A paper on a comparison 
of these theoretical results with experimental measure- 
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ments of the domain wall velocity in CoPt alloy films is 
in preparation. 
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FIG. 1. Hysteresis loop for a lOOA film; T = 300K 



FIG. 2. Domain configuration of a 150 x 150 system during 
the reversal; T = 300K 



FIG. 3. Radius of the domain versus time for the same 
fields as in Fig. 0; T = 0; solid lines are best fitted. 



FIG. 4. Domain wall velocity versus driving field for 
T = 0, 300, and 600K; solid lines are guides to the eye. 



FIG. 5. Scaling plot from Fig. 
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